Non-perturbative electron dynamics in crossed fields 



(N 
O 

o 



in 

(N 



G 

O 

o 



J. M. Villas-Boas, 1 ' 2 Wei Zhang, 1 Sergio E. Ulloa, 1 P. H. Rivera, 2 ' and Nelson Studart 2 

1 Department of Physics and Astronomy, Condensed Matter and Surface Science Program, 
Ohio University, Athens, Ohio 4 5701-2979 
1 Departamento de Fisica, Universidade Federal de Sao Carlos, 13565-905, Sao Carlos, Sao Paulo, Brazil 

(Dated: March 15, 2009) 

Intense AC electric fields on semiconductor structures have been studied in photon-assisted tun- 
neling experiments with magnetic field applied either parallel (-Bn) or perpendicular (B±) to the 
interfaces. We examine here the electron dynamics in a double quantum well when intense AC 
electric fields F, and tilted magnetic fields are applied simultaneously. The problem is treated 
non-perturbatively by a time- dependent Hamiltonian in the effective mass approximation, and using 
a Floquet-Fourier formalism. For By = 0, the quasi-energy spectra show two types of crossings: 
those related to different Landau levels, and those associated to dynamic localization (DL), where 
the electron is confined to one of the wells, despite the non-negligible tunneling between wells. By 
couples parallel and in-plane motions producing anti-crossings in the spectrum. However, since our 
approach is non-perturbative, we are able to explore the entire frequency range. For high frequencies 
u), we reproduce the well known results of perfect DL given by zeroes of a Bessel function. We find 
also that the system exhibits DL at the same values of the field F, even as By 7^ 0, suggesting 
a hidden dynamical symmetry in the system which we identify with different parity operations. 
Symmetries under general parity operations explain many of the features in the spectra, and their 
overall behavior under magnetic field. The return times for the electron at various values of field 
exhibit interesting and complex behavior which is also studied in detail. We find that smaller ui 
shifts the DL points to lower eFd/huj ratios, and more importantly, yields poorer (less effective) 
localization by the field, while other states change also physical character. We analyze the explicit 
time evolution of the system, monitoring the elapsed time to return to a given well for each Landau 
level, and find non-monotonic behavior for decreasing frequencies. 

PACS numbers: 72.20.Ht, 71.70.Di, 73.40.Gk 

Keywords: dynamic localization, tilted magnetic field, double well 
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I. INTRODUCTION 

The dynamics of charged particles in semiconductor 
quantum well structures subject to a time-dependent 
external r-^9|^^f^n$f[W^^H l i s ~| J^^en a subject of intense 
researchBJ3'IBBlal3'iy'yyil3'y The progress of tech- 
niques in nanoscale lithographs and the development of 
the free-electron lasers (FEL)lil which can be continu- 
ously tuned in the terahertz (THz) range, have made 
possible the systematic study of effects only present 
in intense alternating fields in this domain. Examples 
of these effects are the coherent suppression of tunnel- 
ing (or dynamic localizati on) despite interwell tunneling 
in the struct^e,II5ll3ll30 the collapse of minihands in 
superlatticesjl3 absolute negative conductanceJ3'E£l and 
photonpassisted-tunncling (PAT) in resonant tunneling 
diodes£3 the AC Stark Effect ES and many others. These 



effects haw? 
verified O 



^ngrcdicted and partially experimentally 



Intense AC fields on semiconductor structures have 
been also studied in connection with photon-assisted tun- 
neling experiments with magnetic neld-,applied either par- 
allel (-By)] 2 !! or perpendicular (B±)f3 to the interfaces. 
Quantum wells inpti^d-magnetic fields are also of great 
current interestnjCd'ElrL 2 ] as experimental probes of the 
transition to quantum chaos in a mesoscopic system, in 
the presence of a static electric field. In this work, we 
study the combined influence of tilted magnetic field and 



strong AC fields and find quite a rich behavior for differ- 
ent parameters in the problem. 

A double quantum well (DQW) can be treated to low- 
est approximation by considering only the first level in 
each well, resulting in a two level system with a splitting 
A due to the inter-well tunneling. This two-level dynam- 
ics, moreover, plays an important role in understanding 
the behavior of more elaborate dynamical systems, which 
explains why this subject has been so intensely studied 
(for a comprehensive review see Ref. |l3| ). 

In this paper we address the problem of a double quan- 
tum well in simultaneous tilted magnetic field and an in- 
tense AC drive in a non-perturbative approach. To do 
so we make use p£ the Floquet-Fourier formalism devel- 
oped by Shirley£3 For strong oscillating ifclds, the ap- 
proach based on the Floquet theory,E3E3o has proved 
most useful in treating the dynamics of such systems ,-aU 
though some analytical solutions have been proposed. EIEHj 
Shirley's approach is very convenient since a time depen- 
dent problem can be mapped onto a time-independent in- 
finite matrix eigenvalue problem (eventually truncated at 
the desired accuracy) , which yields a series of quasiencr- 
gies describing the time evolution of the system. This 
formalism is also useful as the Floquet states possess ei- 
ther odd or even dynamical symmetry. As a result, to 
explore the effects of magnetic fields, we use a general- 
ized parity operator and its eigenfunctions to provide un- 
derstanding about crossings in the quasienergy spectra. 
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We are specially interested in the dynamics of the system 
in the presence of the intense AC electric field and how 
it changes with the applied tilted magnetic field. Unlike 
most of the previous studies of AC-driven systems, where 
the main attention was paid to the high frequency regime, 
our approach includes the full range of frequency values 
(within the single quantum well level approximation). 

We begin in section II with a complete description of 
the model used, discussing some important properties 
like the parity of the system. We find in fact that this 
system can be understood using two kinds of parity oper- 
ators, one when there is no parallel magnetic field applied 
(S'p) and a full parity operator in the other case (S p ). 
Based in the two level approximation, we also develop an 
analytical description that yields the degeneracy in the 
quasienergy for either the high frequency (A/fiw -C 1) 
or weak electric field (eFd/huj <C 1) regimes, and it is in 
complete agreement with our numerical calculations. 



In section III we make a systematic study of the system 



exploring both the quasienergy spectra and the explicit 
time evolution. We find that in the absence of a parallel 
magnetic field, as the Landau levels are conserved, the 
problem can be understood in terms of two level systems 
for each Landau level. When the parallel magnetic field is 
turned on, however, the crossing levels which originated 
from the Landau level conservation (for Bu = 0), start 
mixing and develop anticrossings. Thus, the symmetry 
changes from S' p — > S p due to B^, lending a transition 
from crossing to anticrossing of levels as the AC field 
amplitude is varied. For all fields, however, we find a 
number of crossings, normally associated to dynamic lo- 
calization, which show no crossing-anticrossing transition 
and little change in dynamics when Bu is turned on. We 
have also performed a complete study of the change in 
the field condition for dynamic localization, as function 
of the ratio A/hio. We find a monotonic drop in the field 
values needed to achieve localization, as the ratio A/fku 
increases (smaller frequency), with successive low field 
DL points disappearing as u> decreases. We are able to 
provide an empirical equation that fits surprisingly well 
the numerical results, although an analytical description 
valid in the different limits is also presented. An analysis 
of the 'quality' of DL for various AC fields is provided by 
a study of the minima in the probability of a particle ini- 
tially in the left well to remain there. This result shows 
a surprising and non-understood behavior. As expected 
for the high frequency limit, the dynamic localization is 
very well defined (fully localizing the particle in a we 
and at fields given by the zeros of the Bessel function) .1 
However, for low frequency, the dynamical localization 
becomes poorer, although the particle never fully leaves 
the well, and presents small 'revivals' for the first and 
third dynamic localization points. This behavior is not 
intuitive and not reported previously. The full range of 
complex behavior is however experimentally relevant, as 



II. THEORY 

The Hamiltonian for an electron in a double quantum 
well structure under a tilted magnetic field and a strong 
AC field can be written as 



H = Hq + eFz cos(a>£), 



(1) 



where Hq = (P — eA) 2 /2m* + Vq is the Hamiltonian for 
the non-driven system. We have chosen the growth di- 
rection along the z-axis and the direction of the magnetic 
field in the x-z plane B = (Bu,0,B±). Using the gauge 
A = (0, B±x — Buz,0), and making use of the transla- 
tional symmetry in the y-direction, the Hamiltonian for 
the non-driven system can be written as 



Hn = 



P 2 , 

2m* 
e 2 B 2 



P 1 

2m* 



e 2 B\\B 



2m* 



nr 



(2) 



where lj±_ = eB±/m* and x' = x+hk y /eBj_, with k y con- 
served throughout. In the absence of Bu, and allowing 
for only one state in each well of the DQW, Eq. (0) is ex- 
actly soluble and its eigenenergies are huj_(n+l/2)±A/2 
(harmonic oscillator in plane with the symmetrical and 
antisymmetrical solutions of the DQW), where A is the 
splitting due to the tunneling. Inclusion of additional lev- 
els in each quantum well is straightforward, and it does 
not add to the qualitative behavior and conclusions de- 
scribed here (except for possible jaulti-photon resonances 
at high field, the AC Stark effectEJ). Moreover, thin wells 
will likely have others levels at much higher energies, and 
our treatment here will be quantitatively accurate. 

The dynamics of the system is governed by the time- 
dependent Schrodinger equation 



(3) 



Since H is periodic in time (H(t) = H(t + r), where 
t = 2tt/uj is the-spried) we can make use of the standard 
Floquet theoryfc-jE3e-3 and write the eigenfunction \tfj) as 



= exp(— iet/h)\u), 
which allows us to rewrite Eq. (|J) as 
{H-ihd t )\u)=e\u), 



(4) 



(5) 



we discuss in section [V 



where \u(t)) = \u(t + r)) is also periodic in time, 
with the same period r, and e is a real-valued param- 
eter termed rite Floquet characteristic exponent, or the 
quasienergy.£3 The term quasienergy reflects the for- 
mal analogy with the quasi-momenturu_fc, characteriz- 
ing Bloch eigenstates in a periodic solidBj Following this 
analogy we can see that |u( m )) = exp(imujt)\u) is also 
solution of Eq. (JsJ) , with eigenvalues e m = e + mftw, if 
m is an integer number m = 0, ±1,±2, .... Then, by 
a subtraction of a suitable integral multiple of hu>, the 
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quasienergy e can be mapped into a first 'Brillouin zone' 
obeying: hw/2 < e < tkv/2. For vanishing intensity of 
the AC field F, the Floquet states are connected with the 
stationary states \a^) of H by \u^) = exp(ilut)\a^) , 
and the quasienergies with the unperturbed states by 
si = E^+lhuj, where the index I counts how many AC en- 
ergy quanta have to be subtracted from the unperturbed 
energy E^ in order jta arrive at the first Brillouin zone 
(the 'photon index') Jiij ._. 

Following the procedure developed by Shirley,Ej we 
solve Eq. (^|) by expanding the function \u) in a Fourier 
series for the temporal component, and linear combina- 
tions of the -B|| = solutions of Hq for the spatial com- 
ponent, 



Cs,n,m\s, Tl, 77l), 



(6) 



where |s) refers to the symmetrical (\S)) and antisym- 
metrical (| A)) eigenfunctions of the DQW, \n) is the Lan- 
dau level function, and |m) is a Fourier expansion 'state' 
((t\m) = exp(imwi)). This leads to the time-independent 
infinite matrix eigenvalue equation 



E 



E s + hcu±(n + -) — e + mhuj 



$s,s' $7i,n' 



e 2 B 2 



l z 2 S 



2m* 



where z 2 s , = (s\z 2 \s'), z s>s , = {s\z\s'), x' nn , = (n\x'\n'), 
and E s = ±A/2. 

A. Floquet states and parity 

The quasienergies in Eq. (p) can be regarded as the 
eigenvalues of a stationary problem analogous to a time- 
independent Schrodinger equation, with \u) playing the 
role of stationary states of the operator H = H — ifodt- 
As a consequence of [H, S p ] = 0, where S p is the parity 
operator defined by S p : (x' — * —x',z — > — z, t — * t + 
t/2X . the Flo quet states |u) have even or odd parity under 
S p Since our expansion basis (Eq. (||)) has also a 

definite parity under S p , the matrix (^) can be separated 
into two blocks for subspaces of different symmetries, one 
for s+n+m = even and another for s+n+m = odd. This 
result is very useful when interpreting the quasienergy 
spectrum, which contains multiple crossings (that come 
from different symmetries), as well as avoided crossings 
associated with mixing between same-parity states. 

It is important to note also that for vanishing Bh , the 
parallel and perpendicular motions are completely decou- 
pled, and in this case the problem reduces to solving a 
DQW system under intense AC field for each in-plane 
Landau level separately. In this £?j_-only case, the parity 



operator can be defined as S' p : (z — > — z, t — > t + r/2), 
and as [H, S' p ] = 0, this yields separate parities under the 
condition s + m — even or s + m = odd. 



B. Evolution operator and degeneracies 

In the Floquet representation, the time-evolution op- 
erator U(t,t ), defined by \tp(t)) = U(t,t )\ip(t )), with 
U{to,to) — I, where / denotes the identity operator, can 
be expressed as 

cxp(—iei(t — t )/h), (8) 

i 

as it is easily checked with the help of Eq. (|B|) . A full pe- 
riod U(t, 0) operator is all we need to construct a discrete 
quantum map, propagating an—initial state over multi- 
ples of the fundamental periodJlj The evolution operator 
formulation allows one to explore the condition of quasi- 
energy degeneracy, which we will see is an important nec- 
essary condition for dynamic localization (although not 
sufficient). 

As mentioned above, without the in-plane magnetic 
field, the problem reduces to a two-level system for each 
Landau level driven by an external time-dependent elec- 
tric field, which has been extensively studied in the lit- 
erature (for a review, see Ref. 13). Although this system 
looks simple, the physics in it is rich, and a complete an- 
alytical solution is either lacking or too-complicated to 
use in the different parameter regimes.u An analytical 
extension to th£-solution in the high frequency limit was 
recently found,c3 providing a perturbative connection to 
this case. Here we give an analytical solution applicable 
in both the high frequency and weak field limits. For 
simplicity, let us consider the system in the absence of a 
magnetic field, so that the Hamiltonian can be written 
in a two level basis as 



H = H 1 +H 2 



A 



eFd 



a x cos(wt), 



(9) 



where a z , a x are the usual Pauli matrices. 

The time evolution operator U satisfies the Schrodinger 
equation in the operator form 



ih^- = HU. 
dt 



(10) 



The difficulty in solving the problem lies in the non- 
commutation between Hi and Hi- For only Hi or 
H2, the solution is Ui = ex-p(—iAa z t/2ti), or U% = 
cxp(— i J Q eFdcos{uit)a x dt/2h) 1 respectively. We write 
the time evolution operator in the form U — UqU' = 
U1U2U', with U' satisfying the equation 



dU' 

ih— = H eS U', 
dt 



(11) 



where H c a = U ^H^U^ — Hi- In the high frequency limit 
(A/huj <C 1) or weak electric field regime (eFd/hui <1), 
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we can write 

U' = cxp (-- J H cS {t')dt^j . (12) 

By diagonalizing U(t), we can obtain the quasienergies. 
After some tedious but straightforward algebra, we find 
that the condition for degeneracy of quasienergy levels, 
(A c s/huj = (mod(27r)) is given by the condition 

A [ 2lr ( A A feFd . . ,A 
— / cos —t' cos — - sin(t') d£ 
rhi) J \nju) J \ rku J 

A /A \ 

+ — 2vr-sin — 2tt =0(mod(2vr)). (13) 
noj \Huj J 

Since the energy splitting A Q g is due to tunneling and 
level mixing by the AC field, this quasienergy degeneracy 
means a suppression of tunneling, or in other words, this 
is a condition for the dynamic localization to occur. It 
is easy to see that in the high frequency limit, Eq. Jl3| ) 
yields the well-known condition J (eFd/huj) = 0. In the 
weak field limit, one obtains the result A/huj = n, n being 
an even integer number (for n odd there is only solution 
for F — and does not lead to dynamic localization). 
As we will see later, these limiting degeneracy conditions 
are in fact realized in our numerical work, and allow one 
to understand the physics of dynamic localization. 

III. RESULTS AND DISCUSSIONS 

In order to perform numerical calculations, we choose 
a symmetric double quantum well (50 — 40 — 50 A 
GaAs/Alo.3Gao.7As/GaAs) with effective electron mass 
m* = 0.067 mo (mo is the electron mass), d = 2(S\z\A) ~ 
90 A is the mean separation between the left and right 
wells, and A = 8.87 meV.cil We have performed a sys- 
tematic study varying the external parameters (Bu, B±, 
u> and F) and analyze the results for different ranges 
of the ratio A/hu>. It is important to emphasize that 
our numerical model is valid for all values of the ratio 
A/hio. However, as will become clear in the discussion 
below, the behavior of the system is markedly different 
for A/hoj < 1 and for A/Huj > 1. We first discuss the 
first regime, A/hio < 1, as it is perhaps more intuitive, 
to later look at lower frequencies. 

A. High frequency, A/hcu < 1 

We start our analysis by considering the quasienergy 
spectrum for two Landau levels and varying Bu . In Fig. 
|l|(a) we show the first Brillouin zone of quasienergy in 
units of photon energy, for A/huj = 0.5 and B± = 1.0 
T. The labels £? SjT j, m for each level, correspond to the z- 
symmetry (S or A), the Landau level index n, and the 
'photon index' m, all at F = 0. For Bn = the par- 
allel and perpendicular motions are completely decou- 
pled and the Landau level number is conserved (only 



transitions with the same Landau level index are al- 
lowed). This explains the level crossings at points such 
as eFd/hoj ~ 1.96 and all others not explicitly labeled, 
according to S' p symmetry (that is, levels crossing have 
different symmetry under S' p ). When we turn on the par- 
allel magnetic field, however, the parallel and in-plane 
motions are coupled, so that the Landau level index is 
no longer a good quantum number. This produces an- 
ticrossings in the quasienergy spectrum at points such 
as eFd/hoj ~ 1.96, while still keeping crossings at the 
points denoted by ?j (localization points), according to 
their S p symmetry. Here we see explicitly that B\\ pro- 
duces a change of symmetry S' p — > S p . Thus the system 
for eFd/hui ~ 1.96 shows a transition from level crossing 
(the states are odd and even in S' p ) to level anticrossing 
(the states are both odd under S p ) due to Bn . It is easy to 
see that there is no such transition for li points (dynamic 
localization points), since at li states have opposite par- 
ity under both S' p and S p symmetry operations. This is 
the essential difference between the point eFd/hu> ~ 1.96 
and li. 

In order to better understand the physical behavior 
of the system at the various crossings, we calculate the 
time evolution of the system for some particular points 
(points indicated by arrows in Fig. |](a)), using Eqs. (Q) 
and (||). Here, the system is prepared in the first Lan- 
dau level of the left well at to, IV'C^o)) = \L,0), where 
\L,0) = \L) ® |0) with \L) = l/V2(\S) - \A)), and |0) 
refers to the lowest Landau level for the in-plane motion. 
We investigate how the probability of finding the particle 
in some state \j3) — \l,n) (I = L,R and n = 0,1,2,...), 
Pp{t) = \((3\ip(t))\ 2 , changes with time. We also investi- 
gate the total probability of the system to remain in the 
left well 

F L (t)=^|(L,n|^))| 2 . (14) 

n 

The probability for staying in the left well at the cross- 
ing/anticrossing point eFd/fuo ~ 1.96 can be seen in Fig. 
|l|(b). We note a large change in its behavior with Bu . For 
= 0, the time evolution shows that the system can be 
found in both left and right wells with the same proba- 
bility and with a rather simple periodic behavior. The 
particle is effectively tunneling back and forth with a pe- 
riod « lOr = 207r/w, and with fast oscillations of period 
r. For Sii =/= the probability presents a more complex 
behavior, with quasiperiodic oscillations. Notice in Fig. 
|l|(b) that the system does not have short-term periods, 
and a long-time plot of Pl (t) in this case yields periods 
w lOOOr and quasi-periodic motion at w 46r. We should 
also point out that the detailed frequency content will 
become even richer if one includes higher Landau levels 
(ignored here for simplicity). 

In Fig. |l|(c) we show results for the point eFd/hu ~ 
2.33, the first localization point l\. We can clearly see 
that in this case, the probability for the system to stay 
in the left well is very high and never goes to zero, in- 
dicating a dynamic localization of the charge, ensued by 
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FIG. 1: (a) First Brillouin zone of quasienergies (in units of huj) as function of eFd/hui (intensity of AC field), for A/hui = 0.5 
and B± — 1.0 T. The label -E s ,n,m refers to the eigenvalues of Ho for vanishing F and B\\ (s indicates the symmetrical or 
antisymmetrical state; n, the Landau level index; and m counts how many photon quanta have to be subtracted or added to 
the unperturbed energy to arrive at the first Brillouin zone). The main difference between Bii = (full line) and Bu = 1.0 T 
(dotted line) are the anticrossings in points such as eFd/hui ~ 1.96. The points h indicate the dynamical localization points, 
and exhibit no change with Bu. (b) Time evolution of the probability Pt(t) of an electron prepared in the first Landau level 
in the left well to be found in the same well for eFd/hu) ~ 1.96 (time in units of the AC field period r = 2n/ui). (c) Same as 
(b), but for eFd/hu; ~ 2.33, the h point. This shows dynamic localization, since the system can be found with at least 70% 
probability in the left well for all time. The inset figures show the time evolution probability for all states of the system, for 
different values of Bii. Note that for Bu — 1.0 T there is a small probability for the system to be found in the first Landau 
level, Pl,i and Pr,i, but Pl is left unchanged, (d) Same as (c) for eFd/hui ~ 5.48 (h)- The dynamic localization is much 
better defined for this point. 



the AC field despite the interwell tunneling allowed in 
the structure. It is important to note that the parallel 
magnetic field has no effect on the position of this dy- 
namic localization point. However, as can be seen in the 
right inset, there is a small probability for the system 
to be found in the second Landau level, n = 1. At this 
point, the system exhibits small transitions between Lan- 
dau levels, but in the same well, so that the localization 
is basically invariant under the presence of the parallel 
magnetic field (although its return period to the n = 
level is now longer, as seen in the right inset in 0(c)). 

For the crossing point at eFd/hui ~ 5.48 (I2), Fig. 
|l](d) shows similar results, except for the fact that the 
time evolution presents better dynamic localization (no- 



tice minima in Pl > 0.92). We in fact find that in gen- 
eral, dynamic localization points are 'better localized' for 
higher F values, resulting in higher minima for Pl and/or 
weaker Pr maxima, although there, is a few exceptional 
examples, as we will see in Fig. wpn 

Another approach to monitor the dynamic localiza- 
tion is to look for the oscillation period in Pl, defined 
naturally as the time for the electron with initial state 
in the left well to go to the right well (with amplitude 
of 100%) and fully back (with amplitude of 100%). In 
Fig. |l|(b), for Sii = 0, it is easy to see that the oscilla- 
tion period is about 10r as mentioned before. Figure ^ 
displays the return periods obtained from direct analy- 
sis of Pl for different values of the AC field amplitude 
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FIG. 2: Oscillation period as function of eFd/hui (intensity 
of AC field), for A/hu) = 0.5, B± = 1.0 T and By = 0. At 
the points U the oscillation period tends to infinity, showing 
the dynamic localization. Solids circles are the results from 
Eq. ([iq), while the full line is obtained from the direct time 
evolution. 



0.15 



0.1 



0.05 



-0.05 



-0.1 



0.5 

s o 



\7¥ 



B n = 0.0T 
B„= LOT 



A/to = 0.1 
B , = 1.0 T 




20 40 60 80 100 20 40 60 80 100 

_l , *h , I , "Ll , 

2 4 6 8 10 

eFd//zW 



FIG. 3: First Brillouin zone of quasienergies (in units of huj) 
as function of eFd/huj (intensity of AC field), for A/hui = 0.1 
and B ± = 1.0 T. Solid (By = 0) and dotted (B n = 1.0 T) 
lines show quasienergy behavior with in-plane magnetic field. 
The inset figures show the time evolution Pl (t) of the system 
initially in the left well for the points indicated by arrows. 



F, for A/huj = 0.5. Notice that the generic period is 
rather short, except for values of the field which produce 
dynamic localization li. Near these values, the return 
period in Pl quickly increases, suggesting indeed a di- 
vergent period right at the k values, and fully consistent 
with dynamic localization behavior. Note that in Fig. 
[l](c) and (d), for example, we cannot see the slow enve- 
lope which yields a finite return time, since the system 
never goes to the right well with more than 50%, and 
all the oscilations remaining correspond to the high fre- 
quency with period r. 

This oscillation period can also be obtained from in- 
spection of the quasienergy spectrum. As Bu = 0, for 
example, only transitions with the same Landau level in- 
dex are allowed, and the oscillation period O p can written 
as 



where £3 refers to the quasienergy labeled £U,o,o and £\ 
refers to Es,o,o in Fig- 0( a ) (quasienergies that evolve 
from the unperturbed state with the same Landau level 
index). The solid circles in Fig. || illustrate the values 
from Eq. (|l5|), in complete agreement with the direct nu- 
merical calculation. We emphasize that as these results 
are rather intuitive, and expected, they provide assur- 
ance that the return period is a good direct measure of 
the degree of localization. As the time evolution of -Pl(£) 
becomes much more complex for different A /Huj ratios 
(or Bu 7^ 0), one can rely in the physical insights it pro- 
vides. 

For the high frequency limit (Huj 3> A), as we dis- 
cussed before, the localization points tend to the zeros of 
the Bessel Function Jo(eFd/Huj) for a double quantum 
well (two level problem). This approximation is quite 



good, even with a tilted magnetic field, as we would ex- 
pect from the different parity of the crossing states under 
S p , and as we can see in Fig. || for A/ Huj = 0.1. The dy- 
namic localization regime here is also much better than 
at A/Huj — 0.5, as shown in the inset (with Pl m in ^ 0.98 
at li, eFd/Huj ~ 2.40, for example). As we will see below, 
the degree (or 'quality') of dynamic localization decreases 
for large A /Huj, even as the quasienergy levels still cross. 

B. Low frequency, A/fko > 1 

In Fig. ^ we show the quasienergy spectrum for 
A I Huj = 2.2. Up to now, all the spectra analyzed could 
have been described without the inclusion of the photon 
index, since it is zero in all cases. In the present case of 
large A/ Huj, however, we need to be more careful, since 
the first Brillouin zone of quasienergies involves multi- 
ple 'photon replicas' and various levels could be involved 
in crossings for different photon index. To understand 
crossing and anticrossing we need to explore the full sym- 
metry of the states under the parity operator S p . 

In Fig. [5] we choose some quasienergy crossing points 
(indicated by arrows in Fig. ^) and carry out the time 
evolution. For the point at eFd/Huj ~ 3.46, the time 
evolution shows that it is a 'normal' point, namely with 
the particle fully jumping to the right well, so that 
the minimum Pl = 0. This behavior suggests that 
the states crossing have different symmetry which the 
tilted magnetic field seems not to affect. On the other 
hand, the crossing point eFd/Huj ~ 4.63 (I2) is simi- 
lar to the dynamic localization points discussed before, 
in the sense that the probability of staying in the left 
well never reaches zero (although here the localization 
is somewhat poor, with minimum Pl — 0.2). More- 
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eFd/to 

FIG. 4: First Brillouin zone of quasienergies (in units of huj) 
as function of eFd/hui (intensity of AC field), for A/hui = 2.2 
and B± — 1.0 T. The h label refers to a localization point and 
it presents no change with Bu . Others crossing/ anticrossing 
points are also indicated (see text). Notice labels E Stn , m in- 
clude different photon index m. 




FIG. 5: Time evolution of the probability Phit) of an electron 
prepared in the first Landau level in the left well, and to be 
found projected in the same well: a) for eFd/hui ~ 3.46; b) 
for eFd/hLu ~ 4.63; c) for eFd/fiw ~ 5.54. Only (b) exhibits 
dynamic localization, although incomplete. 



over, the point eFd/hw ~ 5.54 is also a 'normal' case of 
crossing/anticrossing produced by the inclusion of paral- 
lel magnetic field, so that the probability of staying in 
the left well changes when we include Bu , and presents 
complex time-dependent behavior. For any value of Bu 
however, the minimum Pl — 0, and the particle fully 
tunnels back and forth between the wells. 

It is important to point out that the crossing of levels 
at eFd/hw ~ 3.46, and similar others, is associated with 
their different parity under S p . For example, at the point 
~ 3.46, the levels crossing are those with -Ea,i,-2 an d 
Es,o,+i a t F — 0- These states have then even and odd 
parity under S p , respectively (as given by s + n + rn), 



which remains true for F ^ 0. On the other hand, the 
crossing (anticrossing) at ~ 5.54, and others, correspond 
to Landau levels conservation at B\\ = (or mixing for 
7^ ), and the associated parity under S' p . 
Another region of interest is A/ hw ~ odd integer. In 
that case, different photon replicas are nearly degenerate 
but bring together states with the same symmetry under 
S p (or S' p ). The same parity allows mixing of the levels, 
producing an interesting anticrossing behavior for F ~ 0, 
as shown in Fig. || for A/hw — 0.9 and A/fku = 1.1. 
There we see that in (a), the states labeled s = S (or 
A) have an upward (downward) dispersion for F ~ 
(similar to all cases with A/hui < 1, such as Figs. [j] and 
g). However, for A/hw = 1.1 in Fig. ||(b), the S and A 
states have an apparent exchange of behavior in the F 
field, which is also exhibited for all A/hw > 1 (as in Fig. 
|]). This apparent switching is due to the mixing (and 
level repulsion) seen near F ~ 0, and the fact that the 
labels are always obtained for the F = limit, where 
the AC fieldjia non-existent, and the spatial symmetry is 
well defined.EJ Notice, interestingly, that the related AC 
Stark effect in atoms and superlattice systemsEj occurs 
for A > hw, but this effect is somewhat different. 



C. Incomplete dynamic localization 

Comparing the cases analyzed so far, one promptly no- 
tices that the 'degree' of dynamic location changes with 
the ratio A/hw. In the limit A/hw -C 1, we saw that 
the li values tend to the zeros of the Bessel function, 
and that the probability of staying indefinitely in the left 
well is approximately 1 at all times (full localization). 
For A/hw = 0.5, Fig. [j], the dynamic localization points 
occur for somewhat smaller values of the AC field, and 
the probability of staying fully in the left well is smaller 
(poorer localization). In fact, the particle is able to in- 
creasingly leak out to the right well for lower frequencies, 
even for strong AC field amplitudes. 

A more complete analysis of this behavior is shown in 
Fig. [?]. The full line in |^(a) shows how the dynamic lo- 
calization points change with the ratio A/hw. We have 
numerically computed this function by looking for the 
crossing points in the first quasienergy Brillouin zone 
that produce localization behavior (li). An interesting 
result, as described above, is that the tilted magnetic 
field has no discernible effect on these points, so that 
this result is identically the same as for a simple two level 
problem. Notice that increasing the A/hw ratio not only 
decreases the F values needed for localization, but even- 
tually suppress a given lj, at even values A/hw = 2j, 
where the unperturbed photon replicas would have de- 
generacy points (see discussion after Eq. ((l^)). The odd 
ratios do not yield level crossing and localization, as the 
S p parity symmetry of nearby states is the same (-Es,o,o 
and E/l o _ i, for example, for A ~ hw), and states in fact 
mix (anticross) as we discuss above in reference to Fig. 
^. We should emphasize that the high frequency results 



8 



1 

0. 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0 

-1 
1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 
-1 



— B , 


= 0.0T 


B , 


= 1.0 T 






= 0.0T 




= 1.0 T 




4 6 
eFd/ha 



FIG. 6: Two Brillouin zones of quasienergies (in units of 
hui) as function of eFd/hu (intensity of AC field), for: (a) 
A/huj = 0.9 and (b) A/huj = 1.1, for B± = 1.0 T. Notice, 
despite the similarity of the pictures, that the labels i? Sjn ,m 
are completely different. 



discussed in the literature and yielding dynamic localiza- 
tion at the zeroes of the Bessel functions, agree identically 
with our numerical results. Similarly, the weak field limit 
yields h points at integer values of A/huj, as given by the 
degeneracy condition, Eq. (p"3|). That equation, however, 
does yield an appropriate estimate of the U only in the 
extreme limits of field and frequency. 

Inspection of Fig. ^(a) suggests parabolic curves. In 
fact, in a fully empirical way we use 



Ei AY 



(16) 



where Xj is the j -zero of the Bessel function, Jo(xj) = 
0. This simple fit to the end points describes surprisingly 
well the curves over the whole range, as one can see in 
the dotted lines in Fig. f7|(a). The fit is so good, that it is 
basically indistinguishable for the first localization point 
li, and it only departs noticeably for 1% and I4. 

In addition to the shift in position of the localization 
points lj, and as noticed in the discussion before, we find 
that the value of A/ hut also changes the degree of local- 
ization of the particle. To further study this point, we 
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FIG. 7: (a) Exact crossing points h in the first Brillouin zone 
of quasienergies which yield localization, as function of the 
ratio A/hui. Full lines are from our numerical calculation; 
dotted lines are empirical fits given by Eq. @. (b) Corre- 
sponding minima in the probability to remain in the left well 
for these crossing points. Larger Pl minima indicate better 
dynamical localization, such as for A/hu) < 1. 



have arrived at a quantitative measure of how good the 
localization is, by looking at the minimal probability of 
remaining in the left well during a periodic oscillation of 
-Pt(i). For long times, a large minimal value of Pl illus- 
trates that the localization is 'good'. A small value of 
this quantity indicates that the particle is able to jump 
or tunnel to the right well with a large probability, even 
if the particle eventually returns fully to the original the 
localization is then 'poor'. Figure ^(b) presents the re- 
sulting presents the resulting minimal Pl for different 
dynamical localization points lj . These curves have some 
surprising behavior: (a) Increasing the ratio A/tno pro- 
duces a rapid deterioration of the dynamic localization 
points for all lj. In fact, for A/hw > 5, all localiza- 
tion points shown exhibit such poor localization (mini- 
mal Pl — 0), that the time-dependent oscillations are 
basically non-localized (even though the F field values 
are still somewhat large, Fig. f7](a)). (b) For a range of 
A/huj, the second crossing point I2 is in fact better lo- 
calized than £3, unlike the situation at high frequency. 
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FIG. 8: Minimal probability to remain in the left well for the 
first crossing points l\, for several values of magnetic field. 
The legend indicates the two magnetic field components; the 
first number (left) refers to the intensity of B± and the second 
(right) to Sii in Tesla. Insets show magnification of different 
curve regions, for clarity. 

(c) For li and /3 (but not for I2 and I4), there is a sud- 
den 'revival' of localization, as the 'bumps' in the those 
curves indicate. This 'revival' behavior is not intuitive 
and its source not yet understood. We have attempted 
to link this to the details of the quasienergy spectra (such 
as possible degeneracies or level anticrossings) . We have 
not been able to identify a clear correlation. 

We have also analyzed to what extent the degree of 
localization is affected by other parameters, such as the 
in-plane field, B\\ . Figure || illustrates the minimal Pl 
for different values of the magnetic field. We find that to 
a great extent, the field (both in-plane and perpendicu- 
lar) produces no change in the main drop of the l\ curve. 
However, the height and peak position of the 'revival fea- 
ture' are more noticeably changed by the magnetic field. 
Although we have computed these curves for a wide va- 
riety of total magnetic field amplitudes and/or angles, 
we find that the curve is not drastically modified, and at 
most yields a suppression of this revival feature by a fac- 
tor of nearly two. Notice, however, that the changes ap- 
pear only when both components of the field are present, 



and are larger for the larger total field attempted. Again, 
the quasienergies show no evident feature that one can 
associate with these changes, and its precise elucidation 
remains a challenge for future work. 



IV. CONCLUSIONS 

We have analyzed a double quantum well system in 
the presence of a tilted magnetic field, and a strong os- 
cillating electric field. We find, as other authors, that 
for high frequency, the otherwise tunneling particle is lo- 
calized in one of the two wells at certain values of the 
electric field amplitude. We have however, studied this 
system over the entire regime of frequencies and field am- 
plitudes. We find, in agreement with recent analytical 
work, that the electric field at which localization occurs 
changes monotonically down with lower frequency and 
that this trend continues far away from the perturbative 
regime of A/hus <C 1. We find however, that the degree 
of localization does not change in a monotonic fashion, 
and even shows a certain degree of recovery or revival, 
for lower frequencies. This rather involved behavior with 
the different system parameters is likely to produce a 
number of experimentally observable consequences. For 
example, the varying degree of charge localization for dif- 
ferent frequency values in a given structure would yield 
variations in the generation power of THz radiation. The 
tunability with Bu , on the other hand, would produce en- 
tirely different frequencies in the generation of radiation, 
away from the anticipated driving harmonics. This be- 
havior would then contribute significantly to the already 
rich assortment of different phenomena observed in these 
systems. 
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